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ABSTRACT 

We study the formation of galaxy clusters in the presence of thawing class of scalar 
field dark energy. We consider cases where the scalar field has canonical as well non 
canonical kinetic term in its action. We also consider various forms for the potential of 
the scalar field e.g, linear, quadratic, inverse quadratic, exponential as well as Pseudo- 
Nambu-Goldstone Boson (PNGB) type. Moreover we investigate situtaion where dark 
energy is homogeneous as well as situation where dark energy takes part in the viri- 
alization process. Wc use the Shcth-Torman formalism while calculating the number 
density of galaxy clusters. Our results show that cluster number density for different 
dark energy models have significant deviation from the corresponding value for the 
ACDM case. The deviation is more for higher redshifts. Moreover the tachyon type 
scalar field with linear potential has the highest deviation from the ACDM case. For 
the total cluster number counts, different dark energy models can have substantial 
deviation from ACDM and this deviation is most significant around z ~ 0.5 ~ 1 for 
all the models we considered. 

Key words: Cosmology: Dark Energy, Thawing Model, Halos mass Function, Shcth- 
Torman formalism. 



1 INTRODUCTION 

Over the last decade, the 
Supernovae Typ e Ia(SNIa) 



iRiess et all 120091). Cosmic 



(CMB) (iKomatsu et all 120111 
lations (BAO) (IPercival et al 



ob servational dat a from 
(|Kowalski et alj l20Qg| ; 
Microwave Background 
Ba ryon Acoustic Oscil- 
200S ) and the large scale 



structure surveys ( Cole et alj 2005 ) have confirmed that 
our Universe at present is going through an accelerated 
expanding phase. Till date, there has been a large number 
of proposals to explain such an accelerated expansion. This 
includes the inclusion of an unknown homogeneous matter 
component having a large negative pressure (cosmological 
constant being the simplest example of such fluid), modi- 
fication of gravity of large scale as well as considering the 
backreaction of small scale inhomogeneities in the matter 
distribution. 

Although inclusion of cosmological constant in the en- 
ergy budget of the universe is a minimal way to explain the 
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late time cosmological acceleration and also allowed by all 
cosmological observations, it is plagued by the fine tunin g 
as well as the cosmic coincidence problems jCarrollllioOll ). 
Scalar field models with generic features can alleviate these 
problems and provide an alternative to cosmological con- 
stant. These dynamical scalar field models of dark energy 
are broadly classified into two categories-fast roll and slow 
roll models dubbe d freezing and thawing mo dels. For details, 
see the reference |Caldwell fc LindeJ 12005*) . Among these, 
thawing scalar field models are particularly interesting as 
they can naturally mimic equation of state very close to 
w ~ — 1 which is preferred by all the observational data. 

On the other hand, information about the abundance 
of collapsed structures as a function of mass and red- 
shift is an importan t tool to study the matter distribu- 
tion in the universe (|Evrard et al.1 120021 ). A large number 
of cluster surveys are ongoing or being planned to be setup 
in near future, e.g., PLANCK, eROSITA WFXT which 
would detect a large number of clusters jVikhlinin et al.1 
l2009bl ). Indeed, the mass functions of g alaxy clusters 
have been measured through X-ray surveys (jBorgani et al.l 
l200ll ; iReiprich fc Bohringer 1 12002| ; IVikhlinin et al.1 l2009alb 
via weak and strong lensing studies (|Bartelmann et ahl 



© 0000 RAS 



2 



1 19981 ; iDahldbOOd ; ICorless fc Kindl2009h, u s ing optical sur 



veys, 



, like the SDSS (|Bahcall et al.l 120031 ; IWenl ._ 
well as through Sunayev-Zeldovich effect 1 Tauberl l2005h 
In the last decade several authors have been involved 
in such studies and have found that the dark energy 
not only affects the expansion rate of the background 
and the distance-redshift rela tion but also the growth 
of structure in the universe (IWeinberg fc Kamionkowski 



2010l) 



20031 ; iFrancis et all J2009l; iMortonson 



s, & Ka 
120091; 



Pace et al 



20ld;lKhedekar et all2010bl ; |Khedekar fc Maiumdarll2010bf : 



ManexaJ^booi ^ nce dark energy is expected to 



have an impact on obse rvables such as clust er number counts 
and lensing statistics (|Evrard et al.l [20021 ). Therefore, the 
studies of galaxy clusters would provide a useful tool to 
constrain the model parameters and would help to infer the 
properties of dark energy by discriminating among the dif- 
ferent dark energy models. 

Although the proper way to understand the effect of 
dark energy in the non-linear regime of structure formation 
is through the N-body numerical simulations, one can also 
use a simpler semi-analytical methods which are in reason- 
able agreement with N-body simul ations. One such ex am- 
ple is the spherical collapse model dGunn fc Gottlll972 | ) us - 
ing the Press-Schechter formalism (|Press fc Schechterll 19741 ). 
To overcome the problem of over-prediction of num- 
ber of low(high) mass halos at current epoch for this 
model, an ellipti c al col lapse model was later proposed by 

ISheth fc Tormenl (|l999l ). 

This paper extents the previous work by iDevi fc Sen! 

(|201ll ) by further studying the mo st general analytic form 
of halo mass function introduced bv lSheth fc Tormenl (|l999l ) 
for thawing type scalar field dark energy models with various 
potentials and examine how does the DE affect the abun- 
dance of CDM halos by measuring the galaxy cluster number 
counts. We consider both the canonical and non-canonical 
form of thawing dark energy models. For the completeness 
of the study we check the results for homogeneous as well 
as inhomogeneous dark energy cases. We discuss how much 
they deviate from the conventional ACDM model through 
the linearly extrapolated density contrast 8 c (z) at the red- 
shift of collapse, the number density of halo mass at some 
particular redshifts of collapse and also through the total 
cluster number counts. 

In a recent paper ICampanelli et all (|201ll ) have stud- 
ied the constraints on different dark energy models through 
clust er number counts. They assumed the CPL parametriza- 
tion (|Linder!l2003l ; IChevallier fc Polarskill200 ll) to model dif- 
ferent dark energy mode l s. But it was recently shown by 
iGupta. Maiumdar fc Sen! (|201ll ) that this parametrization 
does not correctly represent all the thawing class of scalar 
field models. Moreover this parametrization is also not suit- 
able for thawing scalar fields with non canonical kinetic 
term. 

The structure of the paper is as follows: in section 2, we 
introduce the thawing dark energy models for both canonical 
and non-canonical kinetic term. We discuss the background 
evolution for such fields considering the difference type of 
potentials e.i. V = </>, V = 2 , V = e? and V = (j>~ 2 and 
also PNGB type. In section 3, we sketch the derivation of the 



equations to calculate the linearly extrapolated density con- 
trast, <5 c (z), at the collapsed redshi ft. Then, we descr i be the 
halos mass function introduced bv lSheth fc Tormenl (|l999l ) 
and calculate the number density and the total cluster num- 
ber counts for our DE models. In section 4, we discuss the 
deviation in the cluster abundance of dark energy models 
we considered from that of ACDM model and finally draw 
conclusion in section 5. 



2 BACKGROUND EVOLUTION 

In what follows, we consider a flat, homogeneous and 
isotropic background universe driven by non-relativistic 
matter and dark energy of thawing type, i.e. + Sl m = 1. 
These thawing type dark energy models are characterized by 
the fact that in the early universe the scalar field is frozen by 
very large Hubble damping and the scalar field starts evolv- 
ing slowly down its potential at the later time. So, the equa- 
tion of state, w(a) = p$/ p$ initially starts with w = — 1 
and slowly departs from it in the later time. We consider 
both ordinary scalar field with canonical kinetic term as well 
as tachyon type scalar field having Born-Infeld type kinec- 
tic term wh ic h are minimal l y coupled to the gravity sector 
(|Sedl2002a! lb1; lGarousi|[200cl ; lKlusonll2000l ). The equation of 
motions for the canonical scalar field and the tachyon field 
are given by 



and 



< + 3H</>(l 



4> + 3H(p + dV/d(p = 
dV/d(j> 



+ 



V 



(i-<n = o 



(2) 



respectively, where dot represents the differentiation w. r. t 
the cosmic time t and the Hubble parameter, H is defined 
as 

tf 2 =(£f = ^(p m + P,). (3) 

Here a(t) is the scale factor, p m = p m o a i is the back- 
ground matter density and p$ = p0o/(i) represents the dark 
energy density with 

»l 



f(a) 



exp 



1 + w(u) 



du 



(4) 



v 



Defining new variables A = — an( i ^ = 

luf^l (^") ' one can f° rm an autonomous system of equa- 
tions involving two observable parameters fi^, and 7 = 
(1 + w) together with the parameter A. For the canonical 
scalar field, it is given as 



7 = -37(2 - 7) + A(2- 7 )\/37^, 

= 3(i-7)n*(i-n*), 

A' = -V3A 2 (r- 1)V7^- 



(5) 
(6) 
(7) 



For the detail derivations see the reference 
(|Scherrer fc Sen! l2008al ). Here the prime denotes the 
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Figure 1. The redshift evolution of the linear density contrast, S c at the redshift of collapse for the canonical scalar field model with 
various potentials for Q m o = 0.25. Left panel: Homogeneous dark energy models. Right panel: Inhomogeneous dark energy models. In all 
panels, the different potentials V{4>) = <j>, <j> 2 , and <j>~ 2 arc represented by different line types along with PNGB model. The ACDM 
case (black solid curve) is plotted for reference. At high redshift, all the models asymptotically approach to Einstcin-de Sitter limit. 



derivative with respect to lna. We solve this system of 
equations numerically for various potentials providing the 
ini tial conditions for 7 an d A as 7 = 0, Aj = 1. As discuss 
in IScherrer fe Sen! (|2008al ) , one can easily see that smaller 
the value of Ai, more the scalar field evolution similar to 
that of cosmological constant A. To check for the maximum 
deviation from the cosmological constant A, we set Ai = 1. 
For Q,f, we choose its initial value in such away so that to get 
required value at present. The difference types of potential 
we considered are V = <f>, V — (j) 2 , V = and V = <f>~ 2 
and they correspond to Y = 0, |, 1 and § respectively. We 
also consi dered the Pseudo-N ambu Goldstone Boson model 
(PNGB) |Frieman et al.lll995t ) which has been characterized 
by the potential: 



V{<t>) = ■ 



*Wf) + 1] 



(8) 



with / = 1. Similar for thawing tachyon field models, the 
autonomous system of equations is given by 



models and can help to infer the properties of dark energy. 
In what follows we calculate the halo mass function for the 
DE models considered in this paper. 



3.1 Spherical collapse 

As our interests lie in finding the halo-mass of CDM mat- 
ter in presence of thawing dark energy, we need to study 
the perturbation of matter inhomogeneity in order to cal- 
culate the linear density contrast at the time of collapse. 
We consider a spherical region of radius r(t) evolving in a 
cosmologically expanding background. The dynamics of this 
spherical region is essentially governed by the Raychaudury 
equation, 



= -4ttG 



, , 1 \ 1 
W{r) + - I P4, c i + -p m d 



(12) 



7' = -67(1 - 7) + 2 X /3^A(1 - 7)* 
S4 = 3fi (l- 7 )(1-^) 

A' = -v/3^A 2 (l-7)*(r-^) 



(9) 
(10) 

(11) 



where p^ c i and p mc i are the density of scalar field and density 
of the matter inside the cluster respectively. It is easier to 
solve the equations after normalizing at the turn around 
point, so we define new variables: 



where F 



v dV_ 



Here also we set the initial condi- 



tions similarly as canonical scalar field case and consider the 
power law potentials men tioned above. For the d e tail calcu- 
lation s see the references |Scherrer fc Se n 2008a bl lAli et al.l 
I2009D . 



a r 
x = — and y — — . 

a t r t 



(13) 



Subscript "t" denotes the turn around time. Now, the 
equations of background evolution and that of perturbation 
reduce to 



x 2 = Ht 2 Sl m ,t[Q m {x)3 



(14) 



3 HALO MASS FUNCTION 



and 



The studies of the mass function of collapsed objects like 
clusters of galaxies can be used to constrain the cosmological 



+ vyl{x,y) 



(15) 
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Figure 2. The rcdshift evolution of the linear density contrast, 8 C at the redshift of collapse for the tachyon scalar field models with 
various potentials for fl m o = 0.25. The different potentials are represented by different line types. Left panel: Homogeneous dark energy 
models. Right panel: fnhomogeneous dark energy models. In all panels, the ACDM case (black solid line) is plotted for reference. At high 
redshift, all the models asymptotically approach to EDS limit. 



where 



U X y) = { I 1 + MKy))] Clustered DE 

[1 + 3w(x)] f(x) Homogeneous DE 



(for I(x,y) see reference (|Basilakos fc Voglisll2007t ))with 
_ P<p,t _ 1 — f^m.t 



(16) 



Pra,t 



(17) 



Here £ represents the matter density contrast at turnaround 
which is defined as 



c 



pmcl,t 
pm,t 



(18) 



The function r(y) is given by r(y) = r t y = £ 1,/3 a t y and 
f2 m (x) by 



fJ m (a;) 



(19) 



1 + ux 3 f(x) 

In addition, the linear density contrast 5 obeys the equa- 
tion: 

8 + 2^8 = AnG Pm 8 = ^Hin mQ a~ 3 5. (20) 

We calculate the linear over density "5 C " at the epoch when 
the spherical region described by equation ()12p collapses to a 
point, t(z c ) — 2t(zt) by solving equation (|20[) . For the initial 



condition 5, = ( Em£L — 1 



we write 



[-] = [c 1/3 !] =(i-w (2i) 



dy x ->o >i/3/ n o \ 



(22) 



Here /3 is the constant term which can be found by solving 
equations |Q3J and |15} after substituting ||STJl and p2]). We 



neglect the higher order terms in a;. For our case the initial 
conditions on 5i are found out as 



^!/3 + £-2/3 
^1/3 + £-2/3 



v £Mm] S± Clustered DE 

/(at) J a t 

^/(l)j Homogeneous DE 



(23) 

As we are dealing with second-order equations, two initial 
values are required, one for the initial over-density 8i and 
other is the initial rate of evolution, S t . S t is generally set 
to S d = 1CP 5 . But we set 8 t = as the result does not 
change considerably, m is set at the matter-radiation equal- 
ity epoch. 



3.2 Number counts 

With the indication from observations that individual galax- 
ies and cluster of galaxies are embedded in extended halos of 
dark matter, the abundance of CDM halos have been stud- 
ied widely. Theoretically, Press and Schetcher were the first 
to describe the abundance of these CDM halos as a func- 
tion of their mass with the assumption that the fraction of 
the volume of the universe that has collapsed into objects of 
mass M at a redshift z follows a Gaussian distribution. The 
comoving number density of clusters which have collapsed 
(i.e., virialized) at certain red-shift z and have masses in the 
range M ~ M + dM can be expressed as: 



dn(M, z) 
dM 



pmo dln<r(M, z) 



A I 



dM 



/(*(*)) 



(24) 



where f(a(z)) is defined as mass function. The standard 
Press-shechter mass function is of the form : 



f(r;PS) 



2 8 c {z) 
n a(M,z 



■ exp 



5l{z) 



2a 2 {M,z) 



(25) 



Although it provides a good general representation of the 
observed distribution of clusters, due to the discrepancy 
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Figure 3. A comparison between the comoving number density of CDM halos for the different thawing scalar field models with that 
of ACDM model by defining a parameter, a in equation 1321 Upper panels: Homogeneous and inhomogeneous dark energy models at 
redshift, z = 2 in the left and right panels respectively. Lower panels: Homogeneous and inhomogeneous dark energy models at redshift, 
z = 0.5 in the left and right panels respectively. In all panels, the various potentials are indicated by different line types. f2 m o = 0.25. 



of over-prediction (under-prediction) of t he number of low 
(high) mass halos at the current epoch, ISheth fc Tormenl 
(| 1999T) introduced an ellipsoidal model of the collapse of per- 
turbations. This Sheth-Torman mass function gives better 
fits to simulated mass function by reducing this discrepancy 
substantially. It has the form: 

(26) 

where there contains three parameters A, a and p which we 
set into A = 0.322, a = 0.707 and p = 0.3 for all the models 
we considered. The Press- Shechter case is recovered for a — 
1 and p = 0. 

The dispersion of the density field on a given comoving 
scale R, containing mass M = 47tp m o-R(M) 3 /3, is given by 

° 2{R) = ^- J" t° 3 P(k)W\kR)f (27) 

where the quantity P(k) is the power spectrum of density 
fluctuations extrapolated to z — according to linear theory 



and W(kR) is the top-hat window function; 

^*>Kw-w) (28) 

D(a) represents the growth function of linear perturba- 
tio n theory and can b e foun d from equation (|20[) (see fig.4 
in ISchaefer fc Kovamal ()2004h ). We normalize the growth 
function such that D(a) = 1 at the present epoch. 

Assuming that the baryon density parameter Qbo <C 
Qcdm,o, the CDM power spectrum can be approximated 
by P(k) oc kT 2 {k) and we use th e transfer function, T(k ) 
introduced by Eisenstein and Hu ()Eisenstein fc Hulll998ah . 

The normalization of the power spectrum is often ex- 
pressed in terms of erg, the rms fluctuation today at a 
scale of 8 h~ Mpc. In our case, using WMAP-7 data 
()Komatsu et al.ll201ll ) we normalise the power spectrum by 
setting ag(ACDM) = 0.80. 

While calculating the halo mass function, we set the 
other cosmological parameters as Q m o = 0.25, Q^o = 
0.75, h = 0.72, Q B0 = 0.0456 and n s = 1. 

We calculate dn/dlogM as a function of M at a partic- 
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Figure 4. Same as figure [3] for the tachyon field dark energy models. The various potentials are indicated by different line types. 



ular z (z — 0.5 and 2) for both the cases of thawing models 
with various potentials and shown in figures ((3| and ((4]). 

The number of clusters in a redshift interval dz, above 
a given minimum(threshold) M = Mmin is obtained from 
dn(M, z)/dM: 



dN ,. r „_ . , dV(z) 

— (M > M min ) = f sk y 



dM ^-A M > z ) (29) 



where f B k y is the fraction of the sky being observed and the 
comoving volume element is given by 



dV a if \ dr 
—z = ^ T {z) Tz 

r(z) is the comoving radial distance out to redshift z: 



•(*) = 



-f 

Ho Jo 



dz 
H(z' 



(30) 



(31) 



where H(z) is the Hubble parameter. For numerical com- 
putation, the upper limit of integration in equation (|29|l is 
replaced by some finite mass value M m ax- The comoving 
volume element is required since the redshift evolution of a 
physical volume in space is model dependent, i.e, depending 
on the form of potentials. 



4 RESULTS AND DISCUSSION 

In this section, we discuss the results for the linear over- 
density contrast, the number density of CDM halos and the 
total clusters number counts for the models we introduced, 
keeping the ACDM model as a reference since the ACDM 
model is currently the simplest model, fitting all available 
observational data despite of its conceptual problems. 

One of our results are shown in Fig[T]where we show the 
linear density contrast, S c , as a function of collapsed redshift 
for both homogeneous (left panel) and inhomogeneous (right 
panel) scalar field models with various potentials. The sim- 
ilar behavior is being shown for the case of tachyon field in 
figure [2] From figure [T] and [2] it can be seen that the linear 
density contrast S c at z = has a significant deviation from 
ACDM case for homogeneous dark energy models. These de- 
viations are comparatively smaller in case of inhomogeneous 
dark energy. This is true for both ordinary scalar field as well 
as tachyon type scalar field. One may expect this type of be- 
havior. For homogeneous case, the equation of state of dark 
energy is greater than w = — 1 (we are not considering Phan- 
tom models). Hence the repulsive effect of dark energy in the 
background evolution is lesser than the ACDM case. This 
results larger linear density contrast for homogeneous dark 
energy. However, for inhomogeneous case, there is an extra 
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Figure 5. The total number counts N as a function of redshift of Z, integrated over mass (from M min = 2* 10 14 M o ) for a survey area 
of 10,000 square degree. All the models are normalised to the same number density of haloes today. Upper panel: The canonical scalar 
field models with homogeneous (left) and inhomogcncous dark energy cases (right). Lower panels: The tachyon scalar field models with 
homogeneous and inhomogcncous dark energy cases in the left and right panels respectively. The different potentials correspond to the 
different line types. The concordance ACDM model (black solid line) is also plotted for comparison. O m o = 0.25. 



repulsive effect inside the cluster due to inhomogeneous dark 
energy. This reduced the linear density contrast for inhomo- 
geneous DE and brings it closer to the ACDM value. At high 
redshift, all the models asymptotically approach to Einstein- 
de Sitter (EDS) limit. Among the potentials we considered, 
the linear potential shows the maximum deviation from the 
ACDM for both the ordinary and tachyon fields as well as 
for homogeneous and inhomogeneou s case. This is cons istent 
with the results earlier obtained bv lDevi fc Senl (|201ll) . 

Next, we investigate a quantity closely related to obser- 
vations, the clusters number density, dn/dlogM , the number 
of the collapsed objects per unit mass per unit volume. Since 
it is only depends on S c and on the growth factor, no appre- 
ciable differences are expected between the models studied. 
In order to see the significant deviation in the cluster num- 
ber density for the thawing models from that of the standard 
ACDM, we define a new parameter, a such that 



10 Q 



dn/dlogM 



(32) 



(dn/dlogM) acdm ' 
This parameter measures the deviation in the clusters 



number density of any thawing models from that of ACDM. 
Larger the value of a more it deviates from ACDM. The 
behavior of a for ordinary scalar field as well as for tachyon 
type field have been shown in figure ([3]) and (j4]) respectively. 
It can be seen from figure ([3]) that a value is higher for higher 
z i.e. the clusters number density for the object collapsing 
earlier have significant deviation from that of the ACDM. 
This is true for both ordinary scalar field and tachyon field 
in homogeneous as well as inhomogeneous DE cases. The 
difference between the homogeneous and inhomogeneous is 
prominent for the object collapsing at the later time (i.e. at 
redshift z — 0.5). Again if we look further, the tachyon field 
models show larger deviations from ACDM at all redshifts 
than the scalar field with canonical kinetic term. Among the 
potentials we considered, the linear potential again shows 
the maximum deviation for the tachyon dark energy models 
for homogeneous as well as inhomogeneous cases. 

An another important quantity that can be derived 
from observations is the total number counts of halos CDM 
above a given mass in a complete survey volume. There are a 
number of cluster surveys ongoing or being planned in near 
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future, e.g., PLANCK, eROSIT A, WFXT which would de- 
tect a large number of clusters (|Vikhlinin et al.ll2009bh . To 
obtain the mass of the clusters detected, one usually utilizes 
a proxy variable which can be the X-ray flux, SZE flux or 
richness of the cluster. The actual quantity chosen would, 
of course, depend on the nature of the survey. The relation 
between these quantities and the mass depends on the de- 
tailed cluster physics. However, for our purpose, it suffices 
to state that the limiting mass M m i n (z) of the survey will be 
essentially determined by the limits on the proxy variable. 

For predicting the number of clusters would be de- 
tected in a typical cluster survey we set the fraction of the 
sky being observed as / s k y = 0.2424, which corresponds 
to a survey area of 10,000 deg 2 , appropriate for surveys 
like eROSITA and WFXT. The value of Af min would de- 
pend on the flux limit of the survey and hence it will be 
a function of z. For simplicity, we choose a constant value 
Mmin = 2 x 10 14 /i _1 Mq, which roughly corresponds to a 
flux limit (in the energy band 0.5-2 keV) of ~ 10 -14 erg 
cm -2 s _1 . We show the calculated value of the total num- 
ber counts of CDM halos as the function of redshift in the 
Figs (0 for all the models we considered. We discuss the 
results below. 

It is clear that there exists significant difference in the 
number counts between different DE models considered. The 
difference is most significant around z ~ 0.5 — 1. To take a 
specific example, we see from the top-left panel that the dif- 
ference in number counts between ACDM and V = <j> homo- 
geneous scalar field model is ~ 2000 at z ~ 0.7. This differ- 
ence is significantly larger than the statistical uncertainties 
(which would be ~ 100 for the type of surveys we are con- 
sidering) , and hence can be used for discriminating between 
different models. The same conclusions can be drawn from 
other panels as well. 

We should mention that discriminating between differ- 
ent DE models would be limited by our ignorance of other 
cosmological parameters like Q, m ,as,,n a - A proper analysis 
of how to constrain the models would involve error estimates 
of the parameters with a combination of different observa- 
tional data sets, e.g., cluster counts, CMBR, BAO etc. This 
is beyond the scope of this paper and would be attempted 
separately in future. 



5 CONCLUSION 

To summarize, we investigate the cluster abundances in cos- 
mological scenario where the universe is dominated by the 
thawing class of scalar dark energy. We consider both the 
ordinary scalar field with canonical kinetic term as well as 
well as tachyon field with DBI form of kinetic energy. More- 
over, we consider a variety of potentials that can give rise 
to suitable cosmological scenario different from concordance 
ACDM model. To study the formation of collapsed struc- 
tures, we consider homogeneous dark energy as well as dark 
energy scenario where dark energy takes part in the virial- 
isation proc ess. We consider the Pre ss-Schetcher formalism 
modified by ISheth fc Tormenl (|l999h to calculate the mass 
function. Subsequently, we show that there exists significant 



difference in the number counts between different dark en- 
ergy models and the concordance ACDM model. Given the 
fact that a large number of cluster surveys are currently 
ongoing as well as a number of future surveys are being 
planned, this can be a smoking gun to distinguish different 
dark energy models from ACDM. Our present work is one 
preliminary step towards that direction. This work can be 
extended to freezing class of scalar field models as well as for 
models with non minimally coupled scalar fields. This will 
be our future aim. 
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